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CALCULATION OF THE LATERAL CONTROL OF SWEPT AND UNSWEPT FLEXIBLE ' 

WINGS OF ARBITRARY STIFFNESS^ 


By Fbankun W . Dirdebich 


SUMMARY 

method similar to that of X AC A Rep. 1000 is presented 
for calculating the e;ffectiL'eness and the reversal speed of lateral- 
control devices on swept and unswept unngs of arbitrary stiff- 
ness. Provision is made for using either stiffness curves and 
root-rotation constants or structural influence coefficients in the 
analysis. Computing forms and an illustrative example are 
included to facilitate calculations by means of the method. 

The effectiveness of conventional aileron configurations and 
the margin against aileron reversal are shown to be relatively 
low for swept wings at all speeds and for all wing plan forms 
at supersonic speeds. 


INTRODUCTION 

Adequate lateral control constitutes one of the more 
significant design requirements for airplanes. The ability 
of the airplane to enter a roll is determined by the control 
potver and is measured by the maximum available rolling 
moment resulting from lateral-control deflection. A meas- 
ure of the degree of lateral maneuverability is the helix 
angle at the \ving tips corresponding to the highest rate of 
roll; the lateral maneuverability depends both on the control 
pover and the damping in roU. 

The control pover and the damping in roll are affected 
by structural fle.xibility. Control deflection ordinarily gives 
rise to aerodynamic loads which tend to deform the wing 
structure in such a way as to reduce the loads on it and thus 
to reduce the control power. If the dynamic pressure of 
the air stream is sufficiently high, the amount of lift which 
results from the structm-al defoimation may be sufficient to 
nullify the effect of the control deflection. The speed and 
dynamic pressure corresponding to this condition are known 
as the aileron reversal speed or reversal dynamic pressure, 
since at a slightly li igb er dynamic pressure a control deflec- 
tion in a given direction would result in a rolling moment 
in a direction opposite to that of the moment on a similar 
rigid wing. 

The present report is concerned with an analysis of these 
problems for swept and unswept wings of arbitrary stiffness. 
The method is based on the analysis of loading of flexible 
wings presented in reference 1 . Since suitable aerodynamic 
influence coefficients are not yet available for antisymmetric 


lift distributions, aerodynamic-induction effects on the lift 
distribution are taken into account only as an over-all 
correction and as a sli^t reduction of the load at the tip, 
as in reference 1. The method is formulated in such a 
manner, however, that" aerodynamic influence coefficients 
may easily be included as soon as they become available. 

The numerical analysis required in any given practical 
case constitutes an extension of the calculations outlined in 
reference 1. The computing forms for the additional calcu- 
lations required for an analysis of lateral-control effective- 
ness or reversal are included in the present report. Their 
use is described in the section entitled “Apphcation of the 
Method,” which may be read without reference to the 
derivation of the method. The material presented in refer- 
ence 1 which is pertiuent to the present analysis is included 
herein. As an example illustrating the method, the lateral- 
control effectiveness and reversal of the wing considered in 
reference 1 are analyzed in this repoit. The reversal speeds 
of several wings derived from this \ving by shifting the 
elastic axis and rotating the wing are calculated to demon- 
strate some general effects of sweep on the aileron reversal 
speed. 

SYMBOLS 
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aspect ratio (b~/S) 

aeroelastic matrices defined by equations (7) 
and (10) 

aiixiliary aeroelastic matrices defined by 
equations (8) and (11) 

aileron-reversal matrix defined, by equation (18) 
section aerodynamic center, measuretl from 
leading edge, faction of chord 
wing span, inches 

wing span less fuselage width, inches 
chord measured parallel to the air stream, 
inches 

average wing chord, inches (S/b ) . 
section lift coefficient (l/gc) * 

section lift-cmwe slope, per radian 
section pitching-moment coefficient referred 
to quarter-chord point 
wing lift-curve slope, per radian 


I Supersedes XACA RM “Calcolattonofilie Effects of Structural Flexibility on Lateral Control of Wings of Arbitrary Plan Form and Stiffness*' by Fnuiklin W. Dfcdecfch^ 
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effective lift-curve slope for twist distribu- 
tions, per radian 

i /c ■ i /Rolling momentN 
rolling-moment coefficient ^ / 

matrix converting torques due to distributed 
loads to torques due to concentrated torques 
matrix converting bending moments due to 
distributed loads to _bending moments due 
to concentrated loads 

section center of pressui-e due to aileron de- 
flection measured from leading edge of 
cliord, fraction of chord 
bending stiffness in planes pei-pendicular to 
the elastic axis, poimd-inches * 
location of elastic axis measured from leading 
edge, fraction of chord 

dimensionless distance along chord from 
reference axis to section aerodjmamic center 
(e-a) 

dimensionless distance along chord from 
reference axis to section center of pressure 
due to aileron deflection {cpi—e) 
factor proportional to the rolling-moment 
coefficient due to aOeron deflection defined 
bj' equation (16) 

torsional stiffness in planes peipendicular to 
the elastic axis, pound-inches * 
integrating matrices for single integration from 
tip to root 

integrating matrices for double integration 
from tip to root 

first rows of matrices [/'] and [//'], respectively 
integrating matrix defined by equation (14) 
dimensionless parameter 

\(EI)r Cr COS*A / 

wing lift-curve-slope ratio 
running air load per unit length perpendicular 
to the plane of symmetry, poimds per inch 
free-stream Mach number 

wing-tip helix angle 
root-t^vist constants 

djmamic pressure, pounds per square inch 
dimensionless dynamic pressure 

{Cl^ q{b'l2)\ c/ cos A 

V mYr 

reduced d3mamic pressure 

dimensionless dynamic pressure 

Cl^ q{b'[2fCr tan A\ . 

{EI)r cos A / 



s 

total wing area including part of wing covered 


by fuselage, square inches 


trace of matrix [4*]" 

t 

running torque due to air load about axes 


perpendicidar to the plane of symmetry, 
inch-pounds per inch 

w 

fuselage width, inches 


distance between the effective root and the 

-- 

innermost complete section of the torsion box 
perpendicular to the elastic axis, inches 

y 

lateral ordinate measured from plane of 


symmetry, inches 


angle of attack due to structural deformatioiij 


radians 

aj 

angle of attack equivalent to unit aileron 


deflection (^1^-^) ^ 

S 

aUeinn deflection measured in planes ])arallel 


to the direction of flight, radians 

e 

moment-arai ratio (fij/ci) 

V 

lateral distance from ■wing mot, inches 

A 

angle of sweepback (measured to the reference 


axis unless specified otherwise), degrees 

I$p] 

influenco-coejBBcient matrix for wing twist in 


planes parallel to the air stream due to 
concentrated unit loads applied at tlio 
reference axis, radians per pound 

I^t] 

influence-coefficient matrix for wing twist in 


planes parallel to the air stream due to con- 
centrated unit torques applied in planes 
parallel to the air stream, radians pc'r inch- 
pound 

Subscripts: 

c/2 

midchord 

D 

divei’gence 

jV) 

flexible wing 

i 

inboard end of aileron 

0 

outboard end of aileron _ 

V 

damping in roll 

R 

reversal 

■ffo 

revei'sal of unswept wing 

r 

at root or effective root 

rw 

rigid wing 

sub 

subsonic 

spr 

supersonic 

w 

wing exclusive of fuselage 

B 

due to (imit) aileron deflection 

Matrix notation: 

{ } 

column matrix 

1 J 

row matrix 

[ ] 

square matrix 

1 1 

diagonal matrix 
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double transpose of a matrix: first about the 
principal diagonal, then about tbe other 
diagonal 
unit matrix 
0 0 0 0 . . “ 

1 0 0 0 . . 

1 0 0 0 . . 

1 0 0 0 . . 

1 0 0 0 .. 


DERIVATION OF THE METHOD 

ASSUMPTIONS 

The assmnptions made in reference 1 are that all deforma- 
tions and angles of attack are small and that the wing 
deformations either are known in the form of structural in- 
fluence coefficients or can be calculated from simple beam 
theory in conjunction with rotations of a flexible root. In 
addition the assumption is made in the present report that 
the angle between the aileron and the wing is constant along 
the span of the aileron and that, in the absence of suitable 
aerodynamic influence coefficients, aerodastic effects due to 
aileron deflection can be calculated on the basis of a modified 
strip theory. 

AIR LOADS 

The lift on a wing section of unit -vvidth paraUd to the 
direction of flight may be expressed in terms of the loading 
coefficient ccijcr as 

( 1 ) 

In the case of a wing at zero angle of attack with ailerons 
deflected the lift may be considered to consist of two parts : 
one due to structural deformation and one due to aileron 
deflection. The loading coefficients for these two parts are 
best treated separately. 

The part of the lift due to the antisymmetric stmctural 
twist can be written in terms of the local values of twist as 

( 2 ) 

in terms of suitable antisymmetric aerodynamic influence 
coefficients [<2o]. Since no such coefficients are available at 
present, modified strip theory may be used, as in reference 1, 
with some saving in labor but at some sacrifice in accuracy. 
With this modified strip theory 


[ V' 

in 

[i.'i = 


where for subsonic speeds the approximate value of effective 
lift-curve slope Cz^ may be obtained from the equation 

^ A cos A 

A+4i cos A. 

and for supersonic speeds, 

^ 4 cos Ac/i 

~ VMo*C03®Ac/2 — 1 


being the lift-curve slope of the section perpendicular to 
the quarter-chord line at a Mach number equal to Ma cos A. 

The lift due to aileron deflection should be calculated by 
a fairly accurate method (see references 2 and 3). This lift 
distribution may be expressed as 



=asSCi,^ 



where the coefficients { ri‘ 1 are the values of C/ calculated 
for a unit effective aileron deflection caS and divided by Cr . 

The combined lift due to twist and aileron deflection is 
then 




( 4 ) 


The torque per unit width of span is the product of the lift 
per unit width and the local moment arm. For the lift 
due to t^vist the moment arm is Cic; for the lift due to aileron 
deflection the moment arm is —etC. (See fig. 1.) Conse- 
quently, the torque may be written as 


where e is the ratio of the moment arms ej/ci. 

THE AEROELASTIC EQUATION 

Method employing stiffness curves. — ^The lifts and 
torques per unit width given in the preceding section can 
be integrated to obtain bending moments and accumulated 
torques about axes perpendicular and parallel to the plane 
of symmetry and, hence, about axes perpendicular and 
parallel to the elastic axis. Integration of these torques 
and moments yields the twists and local dihedrals, .which 
be combined to yield the desired an^e of attack due to 
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^Aerodynamic center 
! £!asftc axis 



.Center of pressure due 



Figube 1.— Definition of geometrical parameters used In the analysis. 


at tke tip is known; otlienvise, niatiices [/'] aiul [//'] ma>' 
be used in supersonic flow also as a fair approximation. 

Method emplo3dng structural influence coefilcients.- -If 
the angle-of-attack changes due to unit concentrated normal 
loads and torques and [^r] have been determined in 
static tests of the actual structure or calculated by a method 
such as that of reference 4, they may be used in an aero- 
elastic analysis in the manner indicated in reference 1 . The 
pertinent aeroclastic equation for the lateral-control i>roblem 
is . . 

{a, — I (9) 

where is the aeroelastic matrix defined in reference 1 as 

= (0“]+ [•!>;.] [O,'] [^]] (10) 

[■3*] is an auxiliary aeroelastic matrix defined by 

= [4>r] [C,'] (0"| [e1-[<I>p] 10,'] (11) 


Structural defoi-mation a,. (See reference 1.) Tlie result- 
ing equation is 

= — og5[4] I (6) 

where k, the ratio of the effective to the actual %ving lift- 
cuiwe slope, is 

_ 4-1-2 cos A 
^ 4-|-4 cos A 


for subsonic flow and has a value of 1 for supersonic flow, 
5* is the dimensionless dynamic pressure, [4] is the aei*o- 
elastic matrLx defined in reference 1 (for subsonic flow) as 


[-4] Q., tan A)[b'] 






{Ei)r\ w, 

El h'j2 ei^Cr cos* A 


b'/2 


[U') [|] 


(7) 


and [4] is an auxiliary aeroclastic matrix defined (for sub- 
sonic flow) by 


{Q.-Q.^ tan A)[1,T 


+(^tan^ vm" 
^{EI)r \_EI 




‘ J L El b'/2 e,^Cr cos* A ‘ ^ ^ ^ I 


( 8 ) 


For supersonic flow, matrices [i] and [//] are used instead 
of matidces [/'] and [//'], provided that the decrease in lift 


(7* is the reduced dynamic pressure defined by 

t /I 


and the load-conversion matrices [(7a'] and [<?/] are defined 
and given in reference 1 . 


SOLUTrON OF THE AEROELASTIC EQUATION 

Calculation of the control effectiveness. — The aero- 
dynamic loading corresponding to a given aileron deflection 
at a given dynamic pressure can be obtained by writing 
equation (6) or equation (9) in the form 


[[ n — [4]] { a, } = — A-gr* a»fi[ 4 ] 



(12) 


Once tlifi right-hand side of equation (12) is evaluated it may 
be regarded as a set of Imovms which in conjunction with 
the coefficients of the matrix [[1]— Kq*[4]] permits a solution 
for the. imknovTis {a,}. The loading corresponding to 
{a,} may then be obtained from equation (4) and hence the 
net rolling moment due to aileron deflection and the result- 
ing wing deformations: 


Rolling moment=2 (13) 

where 

[ir,\ = [H\\+y[I\\ ( 14 ) 

As indicated in reference 1, the solution of the set of simul- 
taneous equations represented by equation (12) may be 
carried out by any conventional method of solving simul- 
taneous equations or by an iteration procedure, the choice 
of method depending primarily on the preference of the 
computer. 

Calculation of the aileron reversal speed. — -The aileron 
reversal speed can be obtained by calculating the rolling 
moment due to aileron deflection, as indicated in (he 
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preceding paragraph, and plotting it against the speed or 
dynamic pressure. The value or values of the speed or 
dynamic pressure for which the rolling moment is zero 
constitute the aileron reversal speeds or dynamic pressures. 
However, a more direct procedure, similar to that used to 
find the divergence speed in reference 1, can be derived as 
follows. 

The term onS ^ equation (12) can be expressed in 

terms of {a,} by means of equations (4) and (13) since at the 
aileron-reversal condition the rolling moment is equal to 
zero. Hence, 


The values of a* and cps are best obtained from experi- 
mental section data at the appropriate Mach numberilio cos A. 
These values, in terms of commonly available quantities, are 

dci 

dS 

da 


and 


cpa=0.25 


dS 

dci 

Us 




or 


OiS= — 


I 




Alultiplying both sides of this equation by 


vields 


Rl/' 


where 




(16) 


Substitution of equation (15) in equation (12) yields 

= (17) 

where the aileron-reversal matrix [J.b] is defined by 


[. 




The value of the dimensionless dynamic pressure nq* at re- 
versal can be found by iterating equation (17) or by expand- 
ing the determinant of the matrix [Ul— k2.*L4.jJI and setting 
the resulting polynomial in nq* equal to 0, as described in 
reference 1. 


APPLICATION OF THE METHOD 

SELECTION OF THE PAEAMETEBS 

The geometric and structural parameters used in the cal- 
culation of the lateral-control effectiveness and related aero- 
elastic properties of a wing are the same as those used in the 
calculations of the aerodynamic loading described in reference 
1. If the root-rotation constants are different for symmetric 
and antisymmetric loadings, those for antisymmetric loadings 
should be used for calculating the [rij matrix used in this 
report. Similarly, the section lift-curve slope, wing lift-curve 
slopes, and local aerodynamic-center values are chosen for 
the ilach number of interest, as described in reference 1, 
except for the calculation” of the aileron reversal speed, as 
discussed in the section entitled “Calculation of the Aileron 
Reversal Speed.” 


Theoretical thin-airfoil values of these pai-ameters are given 
in figure 2 for subsonic and supersonic .speeds. Insuflicient 
information exists at present to permit im accurate correction 
of these data for finite-span effects in aU cases. In a quali- 
tative sense, the section value of ori is known to be a useful 
approximation to the actual value required in the calculations 
of tbis report, except for wings of very low aspect ratio, for 
which the value of as tends to be somewhat larger than the 
section value. 

. The values of cps on a wing tend to be further rearward 
than the values obtained from section data at subsonic 
speeds. The finite-span values may be' estimated by calcu- 
lating lift distributions for a given aileron deflection both by 
unmodified strip theory and by a rational finite-span method, 
such as that of reference 2 or of reference’s,* if the local lift 
given by strip theory is ^sumed to act at the section value of 
eps, and if the difference in the lifts given by the two theories 
is assumed to act at the local aerodraamic center, the chord- 
wise location of the residtant of these t-nn forces may be 
considered to be the three-dimensional value of cpj. On the 
basis of this approximation e»= — at aU points on the wing 
not covered by the aileron. 
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The values of required in this report may also be 

obtained for subsonic speeds by the methods of references 
2 or 3 by calculating spanwise lift distributions in the form 

^ for aj5=l and multiplying the results by — A 

simpler but less accurate way of estimating these values 
consists in using tlie modified strip theory, which is iilso used 
in this report for the calculation of the lift due tq structural 
deformations, .until suitable aerodynamic influence coeffi- 
cients are available. This approximation implies that the 

elements of ^ stations not covered by the 

aileron and 1 for stations covered by the aileron. However, 
in order to take into account the locatbn of the inboard and 
outboard extremities of the aileron with the relatively few 
stations used in the analysis, equivalent values of aj5 have 
to be used. These values, referred to as equivalent 5 values, 
are given in figure 3. They are intended to give a rounded- 

off distribution of '«'liiefi lias approximately the 

same area and the same moment about the root as the 
unmodified strip-theory distribution. The equivalent 5 
values of figure 3 ]>ertain to actual values of oejS equal to 1 ; 

they apply to ailerons which extend from to the tip but 

can be combined to apply to any aileron configuration. 
Several examples are listed in the following table for the six- 
jjoint method, the values of ajS being 1 and the equivalent 

values being read from figure 3(a) as 0.716 for 0.55 
and as 0.293 for ■pj^=0.95: 


Aileron 

ends 

Case 1 

Case 2 

Case 3 

Case 4 

Case 6 



0.^6 

ao6 

a55 

,'q 

0 

. 

6'/2 

1.00 

LOO 

.95 

1.00 


fr'/2 

[ ] 

according to modified strip tlieory 

Case 1 

Case 2 

Cases 

Case 4 

Case 5 

0 

0 

0 

® . 

a 

1 

.20 

0 

0 

0 

JL 

1 

.40 

0 

0 

0 


1 

.60 

.716 

0 

.716 

i ■■ 

■■.284 

.so! 

1 

0 

1 

;i 

0 

.00 

1 

.m 

.707 

1 

0 


Tlic values of case 3 aio obtained from those of cases 1 and 
2 and the values for case 5, from the ones of cases 1 and 4. 


CALCULATION OF THE MATRICES 

A brief introduction to matrix methods is given in the 
appendix of reference 1; a fuller treatment is given in refer- 
ence 5. The numerical constants required in the method of 
this report arc given in reference 1 and in the present report 
for 6-point and 10-point solutions. The elements of the 
matiiccs [7] and [/'] are given in table I and those of the 
matrices [//] and [//'], in table II. The matrix [/]" is 
essentially the double transpose of the matrix [/] and its 
elements are given in table III. If the structural influence- 
coefficient method is used, the required matrices [^i'] and 
[Ct] may be obtained from reference 1 . 

The matrix [A] is calculated as described in reference 1 
and, if desired, by means of the computing forms given in 
reference 1. The steps in the computation may be siun- 
marized as follow's: 


step 

Operation 

® 

•"•L'S'l 



® 

;|J;tan.A[®] 


IVJ 


[®]+[®]+[®] 

® 



[®][®i 

® 

*[®i 

® 

"* * rii'i 

i'/2 BOS’ A ' 


C®I-[®1 

® 

"">ui 

@ 

[@][@] 


tAl=[®]-[@] 


where the notation [(?)], for instance, refers to tlie matrix 
calculated in step 7. For supersonic speeds four additional 
steps are required: 


© 


© 

[®] [®] 

© 


0 
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To proceed with the calculation of lateral control and of 
the aileron reversal speed, the auxiliary aeroelastic and the 
aileron-reversal matrices are then calculated as shown in 
table IV; the numbering of the steps indicated in the upper 
left corner of each block is a continuation of the numbering 
of the steps required to calculate the matrix [A]. The 
auxiliary aeroelastic matrix is obtained in step 15 and the 
aileron-reversal matrix, in step 20. The value of g required 


in step 16 is obtained by poshnultiplying the row [//'oj 



obtained in the step immediately above step 16 by the 


column 


tel 


On the other hand, the square matrix of 


step 18 is obtained by premultiplying the row matrix [@J, 

which is the same row matrix [II'o\ —1 multiplied by — > 

( c ) ^ 

by the column < t • 


A separate set of calculations from step 14 to step 20, 
inclusive, has to be performed for subsonic and supersonic 
speeds; the steps for supemonic speeds can be labeled steps 
14a to 20a to follow the pattern set in reference 1. For a 
1 0-p6int solution, forms similar to those of . table IV can 
easily be drawn up. For the influence-coefficient method a 
set of computing instructions and computing forms can be 
based on equations (10) and (11) in the same way as the 
instructions and forms discussed in this section are based on 
equations (7), (8), and (18). 

Special cases arise when any or all of the values of ei or of 
62 are zero. If only ei, is zero, the value of 6i at some other 
point can. be used as a reference throughout the analysis and 
the parameter g* can be redeflned accordingly. The first 
column of the matrix [@] is calculated in this case by multi- 
plying the first column of the matiix [0] [/'] by the ratio of 
the value ej, to the reference value of ei. Similarly, if some 
other value of «i is zero, say the nth along the. span, e:^ is 
used as a reference but the nth column of [@] is calculated 
by multiplying the nth column of the matrix [0] [!'] by 





, where and c„ are the values of and c at the 


nth station. 

If Cl is zero along the entire span, some of the computing 
instructions given in this report, as well as the ones given in 
reference 1, must be modified somewhat. In table VI (a) 

of refei-ence 1 the matrix L$a)i can be entered in the 
space provided for the matrix. Some of the in- 


structions of table VI (b) of reference 1 and table IV of this 
report are then modified as follows : 


Step 6 


Step 9 


Step 10 



(El)r.W, Qau , „ 

{OJ)r 6.72 tanA ^ J 




Step 11 

As is 

Step 12 

Omit 

Step 13 ’ 

[®] 

Step 14 

cos* A (EI)r 1 
672 (GJ)r tan A 

Step 15 

lA] = [®]-[@] 


AU other instructions are unaffected. 

If e* is zero along the span, table IV of this report may bo 
modified as follows: 

Step 15 [Jn = [@] 

Step 14 may he omitted in tliis case; all other steps in table 
rV are unaffected. 

Similar modifications must also be made if the influence- 
coefficient method is used. 

CALCULATION OF THE AU.ERON REVERSAL SPEED 

The matrix [An] is iterated in table V(a) to calculate the 
critical value of the parameter icg* and hence tlie critical 
speed. The calculation has to be performed once for sub- 
sonic speeds and, if the airplane is to fly at transonic and 
supersonic speeds, once for supersonic speeds. From these 
ciitical values, from the definition of the parametere k and g*, 
and fixim the lift-curve slope the dynamic pressure required 
for aileron reveisal may he calculated and plotted as a 
function of Mach number. If the actual dynamic pressure 
for the altitudes of interest is also plotted on the same chart, 
the lowest intersection of the reveisal with a true-dynamic- 
pressure line will give the revemal Mach number and dynamic 
pressure at the altitude of the true-dynamic-pressure lino. 

The matrices [.Ar] calculated for the special cases men- 
tioned in the preceding section do not all yield the critical 
value of the parameter Kg*. ^Vhen the value of Ci is zero iri 
the root, the critical value of the parameter *.g* based on the 
reference value of Ci will be obtained. If 6i is zero at some 
other point along the span or if 62 is zero along the entire 
span, critical values of the parameter Kg* will be obtained. 
In the case where Ci is zero along the entire span, iteration 
of the matrix [.<4 r] calculated by following the instnictions 
of the preceding section will yield the value of the parameter 
Kg at reversal. 

In some of these special cases, and possibly in other cases 
as w'ell, the iteration procedure may not converge. In 
those cases the lowest value of the parametem {Kg*) a or 
(k2)r is imaginary, so that there is no physical reversal speed 
corresponding to this value, and the wing under considera- 
tion is likely to be safe against reversal (in Uie speed range 
under consideration). If the lowest value of the parameter 
Kg* has the. sign opposite to that of the value of Ct, (or the 
other value of ei used as a reference) or if the critical value of 
Kg has the sign opposite to that of the sweep angle A, the 
reversal dynamic pressure is negative. In tliat case the 
wing also is likely to be safe against reversal, since a negative 
dynamic pressure cannot be obtained at any real speed. 
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However, if the wing is to operate at dynamic pressures 
which correspond to values of kq* or sff much larger than the 
absolute values of the critical values obtained by iteration, 
the nest higher eigenvalues of the aileron-reversal matris 
may have to be calculated by a method such as that given 
in reference 5, page 143. If the nest higher eigenvalue is 
real and of appropriate sign, it defines the critical aileron 
reversal speed. 

Instead of iterating the matrLs to calculate the 

determinant of the matris [Llj— may be espanded 
and equated ^to zero, as noted in reference 1. The result is 
an equation of the type 

- ■ . +Ci(«2*)+1=0 

where n is the order of the inatris — that is, 6 or 10 in the 
case of a 6-poiat or 10-point solution, respectively. Solution 
of this equation yields n values of the lowest real value 
with the appropriate sign is the one that defines the critical 
reversal speed. Instead of actually espanding the determi- 
nant, however, the coefficients G, a . . . On can be. 
obtained in terms of the traces of the powers of the mastris 
[Ab], the trace of a matris being the sum of the elements on 
its principal diagonal, and the ?ith power of the matris: 
[.4 b] being the matris obtained by multiplying [.4 b] by 
itself n— 1 times. If is the trace of [.4.al”, then 

Ci= — Si 


If the same values of xg* are selected as were used in the 
calculation of the aerodynamic loading by the method of 
reference 1, the [[1]— Kg*[.4.]] matris is already available. 
If, in addition. Grout’s method of solving simultaneous 
equations has been used to solve the simffitaneous equations, 
part of the auxiliary matrix is also available so that calcu- 
lation of the a, values for the aileron loading requires very 
little time. However, the iterative solution does not have 
this advantage. 

In some of the special cases discussed in the preceding 
sections, care must be taken to use the proper parameters in 
conjunction with the matrices calculated for these special 
eases. In the case where Ci,. is zero, the values of Kg* must 
be based on the reference value of ei selected in calculating 
the matris; in the case where Ci is zero along the entire span, 
the parameter Kg must be used instead of Kg* in table V (b). 

The resulting a, values may be added algebraically to the 


values of 



multiplied by ^ and 
^ 0 



as indicated in 


steps 4, 5, and 6 of table V(b), and plotted over the span 

( cc 

to yield the net aerodynamic load distribution 

which pertains to a unit value of agS on the flexible wing. 
The roUing-moment coefficient due to thU forcing loading 
(over both wings) may be obtained from a dimensionless 
form of equation (13) 



C2 — — ^ CGiSi-[-Ss) 

^3== (CzSi-t- CiSj+Ss) 


Gs = — i(G»_iSi-[-C'n_2S2+ - • - +ClSa-i + Sa) 

tlf 

Unless certain types of automatic computing machinery 
are available, the iteration procedure is generally preferable 
to the procedure based on the expansion of the determinant. 


CALCULATION OF CONTROL POWER AND itANEOTERABnJTX 


The calculation of the twist distribution for a given 
aileron deflection may be carried out in table V (b) , which is 
similarto table^TL (b) of reference 1. The matrix [[1] — Kg*[4,]] 


is entered at the left, and the column 


is enter- 


ed at the light. This column is premultiplied by the 
[4l] matrix obtained in step 15 or step 15a and is entered in 
the second column at the ri^t, which in turn is multiphed 
by —Kg* to 3 deld the third column. The simultaneous 
equations with the coefficiaits at the left and the knowns at 
the right (the third column) are then solved for the unknown 
ctf values. If preferred, an iterative solution of the type 
discussed in reference 1 may be used instead of Grout’s 
method (reference 6) for which table Y is set up . A computing 
form similar to that of table "VTI (c) of reference 1 may be 
used for this purpose. 



(19) 


This coefficient, which is a direct measure of the rolling 
power, is seen to be dependent only on g/g® (except for the 

factor' ^ since and (Kg*)/, is constant for a 

given speed range. 

The rolling maneuverability depends not only on the 
rolling power but also on the damping in roll. The rate of 
roll per unit aileron deflection (measured in a plane parallel 
to the plane of symmetry) is given by 



where Cn is the forcing coefficient calculated from equation 
-(19) (if the contributions of the pressures on the fuselage are 
n^lected) and C^, is the damping coefficient calculated from 


equation (19) with a column of values of 



calculated 


by the method of reference 1 for a case where 
modified strip theory is used, the desired column is 





where a, is the structural deformation associated with the 
given values of 
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ILLUSTRATIVE EXAMPLE 

The method described iu the preceding sections has been 
used to analyze the lateral maneuverability of the wing con- 
sidered in the illustrative example of reference 1. The re- 
quired additional parameters of this wing are presented in 
table VI (a), which follows the form of table IV (a). Modi- 
fied strip theory has been used for calculating > • The 

( 

equivalent value of 5 at the station 0.4 is obtained 

from figure 3 for the given values of and The 

auxiliary aerpelastic matrix for the subsonic case has been 
calculated by following thn form of table IV (b ) ; the re- 
sulting matrix is shown in table VI (b). 

The aileron-reversal matrix for the subsonic case is calcu- 
lated by means of the form of table IV (c). Several of the 
steps, as well as the result, are shown in table VI (c) for the 
subsonic case. In these calculations the contribution of the 
matrix [/'ij to the matrix [//'oj has been neglected, so that 
the matrix [IIW has been used instead of the matrix [II' q\, 
a procedure which is not recommended in general. Iteration 
of the aileron-reversal matrix (by means of the form of table 
V(a) or otherwise) yields a value of (Ag*)s=2.364. A 
similar calculation for supcmonic speeds yields a value of 
(<cg*)H=0.1280. From these tw'O values and the definition 
of the parameters k and q* the dynamic pressure required for 
reversal baa been calculated and is plotted against Mach 
number in figure 4. Also shown in figure 4 for comparison 
are the dynamic pressures required for divergence as well as 
the actual dynamic pressures at sea level and at an altitude 
of 25,000 feet. Where the dynamic pressure required for 
reversal is less than the actual dynamic pressure, the aileron 
control is reveised. For the example wing, reversal is likely 
to occur at a Mach number of appi-oximately 1.3 at sea level. 

The aerodynamic loading due to aileron deflection has 
been calculated by means of the form of table V (b). For the 
subsonic case and for k 2*=0.552 the [[I]— Kg*[A]] matrix is 
that shown in table X(b) of reference 1. The three columns 
to be entered at the right of table V(b), as well as the four 
columns obtained at the bottom of table V(b), are as follows; 



1 

1 

® 

|4;1, 

ra{®i 


0 

O'; 

0 

0 

.1^ 

“=1043 

.206 

.6J^00 

.-.2816 

1 

,saii7 

—.-4530 

1 

1.0l(7S 

;-.60n 

1 

i.ote 

-.6882 





FIral 

matrix 

|t®t 

W'®' 

i®t+K) 

!».! 

0 

0 

0 

0 

-.080 

-.003 

-.070 

—.0699 

.067 

.063 

.005 

1999 

.500 

.478 

.001 

-.3392 

.435 

.344 

.662 

-.4477 

.37V 

.300 

.523 

-.4773 

U-f'dl® 1-0.179 




Figuee 4.— EBcctofMucb number on the divergence and alleron^^!vcrsal dynamic pressure 
of tbe eiampic wing. 


The lift distribution due to aileron deflection 



obtained by using the modified strip-theory values of < jr — / 

■ . ( 
is plotted in figure 5. For the flexible wing the lift distribu- 
tions due to the calculated twist distributions, sucli as the 
one shown in the next to the last column of tlie foregoing 
tabulation, must be added algebraically to tlie lift distribu- 
tion due to ailei-on deflection. This addition is best per- 
formed by first plotting the lift distributions due to twist 
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Case 

lliD 

tlQIt 

Condition. 

1 

0 

0 

Subsonic 

2 

-0.26 

0.234 

Subsonic 

3 

-0.60 

0.467 

Subsonic 

4 

-0.C0 

0.677 

SnpKsonlc 

9 

' 

-LOO 

L154 

Eupetsonic 


FinuBE 5.— Loading oC exam^ wfng with. aflerozL deflected . 


separately and then adding them point for point to the lift 
distribution due to aileron deflection. The net distributions 
obtained in this manner for several cases are shovm in figure 

a. Thedistributionforcase 5 ^supersonic speeds, — 1-00^ 

indicates that the wing is operating at a speed above 
its reversal speed; from the given values of s*z> and q*jt the 

dynamic pressure for the case of —=—1.00 can be seen to 

2d 

exceed by 15.4 percent the dynamic pressure required for 
aileron reversal. 

The roUing-moment coefficient is obtained from equation 
(19) or by adding the moments corresponding to the aileron- 
distribution curve and the twist curve algebraically. (As 
stated previously, the contribution of [7'iJ to [J/'o] has 
been neglected in these calculations.) The ratio of the 
flexible-wing rolling-moment coefficient obtained in this 
manner to the corresponding rigid-wing rolling-moment 
coefficient is plotted in figure 6 (a) against the ratio — g/gi>; 
for the rigid wing the value of at subsonic speeds is 

and at supersonic speeds is 

(f7O.„=0.026C.^^ 

The lateral maneuverability is calculated by means of equa- 
tion (20) with the damping coefficients calculated in reference 



(a) VadationwiUidjiiamtcpressareiatla. 

Cb] Variation, with Mach munbeiasdalUtiiiie. 
Fiouaa 6.— Control power and edecttroiess of example wing. 


1 and is also plotted in figure 6 (a) as a fraction of the rigid- 
wing value, which at subsonic speeds is 

and at supersonic speeds is 

Both the maneuverability and the control power become 

zero at a value of -^=—0.87, which is the ratio of the re- 

versal to the divergence dynamic pressure at supersonic 
speeds, as is shown in figure 4. 

Since the ratio g/go has been determined as a function of 
altitude and ilach number in figure 4, the parameters of 
figure 6 (a) can be plotted as functions of altitude and ilach 
number, as has been done in figure 6 (b) . The maneuverabil- 
ity and, to a lesser extent, the control power are relatively low 
at supersonic speeds, particularly at low altitudes. Since at 
high speeds even a small value of p6/2T’’ implies a fairly large 
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value of the rate of roll p, this situation is not necessarily 
alarming. The wing in question should have adequate 
control at ah speeds for altitudes greater than about 20,000 
feet. 

DISCUSSION 

The method of this report is based in essence on a numer- 
ical integration by means of matrices of the differential 
equations of structural equilibrium. The actual stiffness 
distributions, root rotations, and the lift and pitching- 
moment distributions of the undeformed wing can be taken 
into account as accurately as they are known. The commonly 
made simplification of treating the wing as an aggregate of 
constant-chord segments with all flexibility concentrated at 
the ends and aU forces at the midpoint of the segments 
is not resorted to in this report. No time-consuming 
graphical integrations nor trial and error procedures are 
used. The aileron reversal speed is calculated by means of 
an iteration, but each cycle of this iteration consists of a 
single matrLx multiplication so that the entire procedure is 
straightforward in apphcation, and usually the results 
converge rapidly. If prefeired, the iteration procedure can 
be replaced by an expansion of the determinant of the 
matrix [111— as outlined in this report. 

The purpose of this section is to discuss the assumptions 
and limitations of the method of this report and to indicate 
the effect of certain design variables on the aUeinn revei'sal 
speed by means of the results of a few calculations for the 
e.xample wing and some others related to it. 

'assumptions and limitations of the method 

The discussion of the aerodynamic and structural assump- 
tions in reference 1 is also pertinent to the analysis of this 
report. This discussion may be summarized as follows: All 
angles of attack, structural deformations, and control 
deflections must be sufficiently small to give rise only to 
linear aerodynamic and structural forces. When structural 
influence coefficients are used, no further assumptions are 
necessary concerning the structural deformations; when the 
stiffness curves are used, elementary beam theory as cor- 
rected by root rotations must be applicable. If elementary 
beam theory is inapplicable — that is, if shear deformations, 
shear lag, and bending-torsion interaction cannot be neg- 
lected — a more refined method than elementary beam theory, 
such as the method of reference 4, can be used to c^ilculate 
structural influence coefficients whick.can be used in the 
method of the present report. When suitable aerodynamic 
influence coefficients are available, no further assumptions 
need be made concerning the aerodynamic forces; if no such 
coefficients are available, the assumption must be made that 
modified strip theoiy is sufficiently accurate to calculate the 
aeroelastic effects of intei’est. . 

In the present report additional aerodynamic assumptions 
must be made, primarily, because although accurate aero- 
dynamic information can be used in the method of this report 
such information is not available in many instances. For 
instance, no suitable aerodynamic influence coefficients are 
available as yet for antisymmetric lift distributions so that 
modified strip theoiy has to be used for the lift due to 


structural deformation. For the lift due to aileron deflection, 


which is used in the form of the coefficients 



the best 


available information should be used; for unswept iviiigs of 
moderate or high aspect ratio the method of reference 2 gives 
accurate results, and for swept wings or \vings of low aspect 
ratio tlifi...raethod of reference 3, ivith certain modifications 
explained m reference 6, may be used to calculate the desired 
coefficients. However, information concerning the param- 
eters ctj and cpi for wmgs of finite span is very meager; the 
suggest^ means of estimating them give results which must 
be used’wdth some caution. If experimental results arc 
available for these parameters, they can, of course, bo used in 
this method. 

Modified strip theory should not in general be used to 


calculate the coefficients ( 77^ i ' as W'as done in the illnstra- 

tiye e.xample. If it is desired to use this approximation, the 
equivalent 5 values of figure 3 may be used to obtain a 
suitable fairing of the lift-distribution curve at tfie aileron 
ends. However, these equivalent values are premised on the 
use of the [/'] and [//'] matrices and should not be emiiloyed 
for anj' other purpose than that indicated herein. 

Two additional structural assumptions are also made in 
this report. In the first place, the angle b between the wing 
and the aileron is assumed to be constant along the span. 
This assumption appears to have been made in almost all of 
the published investigations into the problem of lateral- 
control reversal and appears to have yielded satisfactory 
results; the shorter the aileron and the greater the number 
of points at which the aileron is supported and at which its 
hinge moment is taken out, the more nearly true is the 
assumption. Also, the control linkage is assumed to be, 
stiff so Jhat the aileron angle for a given stick displacement 
is independent of the dynamic pressure. However, in order 
to accoimt for the control-linkage deflection, it is necessary 
only to calculate the ratio of the true aileron angle at a given 
dynamic pressure to that at zero dynamic pressure for the 
same stick position. The calculated control moment and 
maneuverability must then bo reduced by this factor to get 
values of these quantities for a given slick displacement. 
Since deformations of the control linkage only affect the 
aileron effectiveness, they have no bearing on the revereal 
speed. On the other hand, these deformations may lead to 
aileron divergence for wings with heavily overbalanced 
ailerons. This problem, as well as the problem of wing- 
aileron divergence, has not been considered in the present 
analysis. 

The fuselage and tail do not contribute any appreciable 
amounts to either the control or the damping moment so 
that their effects may ordinarily be neglected for the jmrposc 
of lateralrcontrol calcidations. Similarly, the effect of wing 
camber does not enter into the problem, because the only 
impoi-tant effect of camber is to give the flexible wing a sym- 
metrical lift distribution if it is set at the angle of attack 
which would give zero lift for the rigid wing; this symmetrical 
lift distribution has no effect on the lateral-control problem. 
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As in reference 1, the effects of the iaertia loading on the 
aerodynamic loading hare not been considered explicitly in 
the analysis of this report. As pointed ont in reference 1, 
however, the structural deformations due to the inertia load- 
ing may be calculated conveniently by means of the integrat- 
ing matrices and then be considered as part of the geometrical 
angles of attack. This procedure may be applied in the 
case of a rolling wing to determine the change in rolling 
moment for a unit rolling acceleration at any g^ven Mach 
number and dynamic pressure. This rolling moment must 
be taken into account in estimating the rolling accelerations 
due to a given forcing moment at any time before the steady- 
roll condition is reached. 

At transonic speeds there is considerable uncertainty 
in the aerodynamic parameters. The control power is 
directly proportional to the value of the parameter 
dc 

which may be quite low in the transonic region 


due to the fact that the aileron is located in a region of sepa- 
rated flow. The method of this report is appKcable to this 
case if the value of Cjj is known for the rigid wing and if the 
aerodynamic forces due to aileron deflection and due to 
twist can be superimposed linearly. If, for instance, the 
decrease in this parameter due to flow separation is 40 per- 
cent at a given Mach number and if the loss in control power 
due to wing flexibility amounts to 20 percent, then the total 
loss is 52 percent. However, the loss in maneuverability 
due to the decrease in may be much less than the loss in 
control power, since a decrease in Cjj is usually accompanied 
by a decrease in C; and Ct, and hence in the coefficient of 

damping in roll. ~ 

Should the value of the parameter Cj, decrease to zero or 
reverse, aileron reversal will occur. This type of reversal is 
altogether different from the type of reversal discussed in 
this report since it is due entirely to aerodynamic action, 
whereas the reversal of concern in this report is due to aero- 
elastic action. Both types of reversal are largely independent 
of each other; aerodynamic reversal is likely to occur at a 
given Mach number regardless of the stiffness of the wing, 
whereas aeroelastic reversal wfll occur ordinarily at a 
different speed which is unaffected by the aerodynamic 
effectiveness. 


EFFECTS OF SOME DESIGN VAHIABLES ON THE AILEBON RF. VERS AT. 

SPEED 

Some general effects of sweep and of the moment arms fi 
and €i on the aeroelastic reversal speed may be of interest. 
The ratio of the reversal parameter of a given wing 

to that of the unswept wing obtained by rotating the given 
wing is shown in flgure 7 (a) plotted against a func- 

tion of the sweep angle for subsonic and supersonic speeds; 
the two curves were obtained by considering the wing of the 
illustrative example to be rotated in such a manner as to 


keep the parameters 


Cr cos^A 


and 




as well as the 


672 {OJ)r 

chord, stiffness, and moment-arm distributions and 
constant. 




(a) Effect of sweep on ictks^ paismeter CcO <• 

Cb) Effect ot momerLt-nrm ratio eon. reversal parameter (is*) «tftirmiswcpt wings. 
Ficcanr.— Effects of sweep and moment^rm ratio on reversal speed. 


Both sweepback and sweepforward apparently tend to 
decrease the reversal parameter and hence the reversal speed. 
At supersonic speeds or, more specifically, at small values of 
cos*A ^ {EI)r 


the parameter 


the reversal speed for the 


672 A' {GJ)r 

sweptforward wing is somewhat lower than that of the 
sweptback wing; whereas at higher values of the parameter 
the varia tion o f the reverse speed with the sweep parameter 

tan A is more nearly symmetrical with respect to 

the zero-sweep case. There are some indications that this be- 
havior is not typical of all wings but rather is due to the fairly 
large variation of the values of ei, e^, and e over the span 
of the example wing. In general it appears that, for sTna.n 

^2 cos^A. l( ?? 

values of the moment-arm parameter '' ^ ' 

the variation of the reversal speed with the sweep parameter 

[WJ)r 

y(jEi)r 

large values of the moment-arm parameter, the reversal 
speed should tend to be lower for sweptback wings than for 
sweptforward wings. 

The variation of the reversal speed of an unswept wmg 
with the moment-arm ratio is shown in figure 7 (b) for wings 
which have the same distributions of the parameters Ci/fi, 
and along the span but have different values of ej ^ and 


tan A 


should be nearly symmetrical and that, for 
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The parameter (k 2 *)j!o is plotted against the ratio 

where the value of e is selected at the midaileron station. 
It is seen that the plot is linear for both the subsonic and 
supersonic case. The difference in these cases is due to 
the different variations of and e* along the span; if the 
variations were the same or if ei and Cj were constant along 
the span, the two Hues of figure 7 (b) would coincide. Since 

the reversal parameter ^ proportional to and 

since the reversal djuiamic pressure is directly proportional 
to the reversal parameter and inversely proportional to the 
value (by definition of the parameter nq*), it follows 
that the reversal dynamic pressure is approximately propor- 
tional to the ratio ^ From figure 1 it is seen that 

the sum of Ci and represents the distance from the aero- 
dynamic center to the center of pressure of the lift due to 
aileron deflection and is independent of the location of the 
elastic axis. This fact corroborates the commonly made 
observation that the reversal speed is independent of the 
location of the elastic axis in the case of unswept wdngs. 

The control power and maneuverability cannot be related 
to the structural and geometric parameters in as relatively 
simple a manner as the reversal speed. The control power 
is a fimction of both the ratio qlqn and the ratio Wgij; it 
normally decreases with qlq^, the rate of decrease being 
slow at first and then more rapid for positive values of 
2 a/?i 3 (which would generally be obtained for unswept and 
sweptforward wings) and being rapid at first and then slower 
for negative values oi qslio (which would generally be 
obtained for sweptback wings). The variation of the 
maneuverability should generally be similar to that of 
the control power since the damping coefficient decreases 
(or in the case of unswept and sweptforw'ard wings increases) 
steadily with q/q^ and is independent of gs/'an- 

From the calculations for the example wings it appears 
that the control power and maneuverabUity of sweptback 
wings tend to be relatively low, particularly at supersonic 
spe^s. A combination of high sweep and large moment 
arm Ca may lead to an undesirably bw maneuverability. 
Of course, any increase in the purely aerodynamic effec- 
tiveness at of the aileron-airfoil combmation results m a 
proportional increase in the lateral-control effectiveness. 
At supersonic speeds as is proportional to the aileron-chord 
ratio cjc, so that an increased aileron chord results m 
greater maneiiverability ; at subsonic speeds an increase 
in the aileron chord is less effective. Another obvious 
means of raising the reversal speed and of mcreasing the 
control power is to increase either the torsional stiffness 
or the bending stiffness of the structure. In some cases, 
however, the increase of the reversal parameter («?*)« due 


to a cliange in the parameter tan A (see fig- 7 (a)) 

produced by a decrease in the torsional stiffness (GJ) r may 


be so rapid as to cause a net increase in the reversal speed. 

If the sweep, the moment arm Cj, the stiffness, and the 
aileron effectiveness cannot be changed sufficiently to 
mcrease the maneuverability, it may be necessarj' to resort 
to uncoventional control devices. Leading-edge ailerons, 
for instance, have negative values of the moment arm e*, 
so that wings equipped with them tend to reverse at verj 
high speeds, if at all. This type of configuration has the 
additional advantage of relatively high effectiveness at 
transonic speeds. The effectiveness of leading-edge ailerons 
at subsonic speeds is so low', however, that they would have 
to be used in conjunction with trailing-edgc ailerons to 
assure satisfactory lateral control at low' sul;sonic speeds; 
furthermore, they pose some other aerodynamic as well as 
structural and mechanical problems. Similarly, from an 
aei-oelastic point of view spoilers appear attractive because 
they tend to have small or negative values of Cj, but they 
also pose certain design problems. Consequently, both 
these devices require careful consideration before they are 
used to alleviate aeroelastic difficulties in any specific case. 

CONCLUDING REMARKS 

A method has been presented for calculating the effective- 
ness and the speed of reversal of lateral control as well as of 
the aerodynamic loading and the rolling moment jirotluced 
by aileron deflection on swept fle.xible wings of arbitrary 
stiffness. 

It has been shown that the aileron reversal speed decreases 
with both sweptback and sw'eptforward wings and that the 
effectiveness of conventional aileron configurations on sw'cpt- 
back wings at supersonic speeds tends to be relatively low. 
The control. effectiveness and the resulting maneuverability 
of the airplane may be increased by varj'ing some of the 
design parameters such as the structural stiffness and, if 
necessaiy, resorting to unconventional control devices, such 
as leading-edge ailerons or spoilers. 

Langley Aeeonadtical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., AprU 6, 1951. 
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TABLE I.— VALUES OF THE INTEGRATING MATRICES [I] AND [/'] 

(a) Six-Point Solution 



(b) Ten-Point Solution 
[/'] 
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TABLE II.— VALUES OF THE INTEGRATING ^MATRICES [//] AND [//'] 

(a) Six-Point Soldtion. 

[//] un 


n 

6'/!2 

0 

.3 

.4 

.0 

.8 

.0 

0 : 

0 

a05333 

0.05333 

■ 0.16000 

0.073M 

0.13702 

. 2 .. 

-.00107 

.01000 

.02500 

.1OG0T 

.05418 

.1077S 

.4; 

0 

0 

0 

.05333 

.03581 

.07758 

.6 

0 

o‘ 

-.00167 

.01000 

.01548 

. 04741 

,s 

0 

0 

0 

0 

-.00152 

.01724 

.g 

0 

0 

0 

0 

-.00108 

.00410 


1J 

b '{2 

1 

.2 

.4 

.6 

.8 

.0 

1.0 

0 

1 

0 ! 

1 ). 06333. 

0.06333 

0. 16003 

0.08000 

ai2ooo 

0.03333 

.2 


ooloj 

..01060 

.02500 

.10607 

.'00000 

.06333 

.02007 

.4 

0 ; 

0 

0 

.05S33_ 

.04000 

06667' 

.02000 

.0 

0 : 

0 

-.00167 

.01000 

.01833 

.04000 

.01333 

.8 

® ! 

0 

0 

0 

0 

.01333 

.00067 

.g 

0 

0 

0 

0 

-.00042 

.00250 

.00262 

1.0 

0 

1 H 

0 

0 

0 

0 

0 

0 


(b) Ten-Point Solution 


1? 

6'/2 

0 

•1 - 

.2 

.3 


•6 V 

.6 

.7 

A 

.g 

0 

0 

0.013333 

0.013333 

a 040000 

0.020667 

0.066667 

0.040000 

0.093333 

0.046476 

0.137920 

.1 

-.000417 

.002500 

.006251 

.026667 

.020000 

.053333 

.093333 

.080006 

.040476 

.122&U 

.2 

0 

0 

0 

.018333 

.013333 

.040000 

.026007 

.066067 

.034477 

.107750 

.3 

0 

0 

-.000417 

.002500 

.006251 

.026667 

.020000 

.053333 

.028476 

.092666 

.4 

0 

0 

0 

0 

0 

.013331 

.013333 

.040000 

.022476 

.077580 

.5 

0 

0 ., 

0 

0 

-.000417 

.002500' 

.006251 

.020007 

.010477 

.062405 

.6 

0 

0 

0 

0 

0 

0 ll 

0 

.013333 

. 010476 

.047410 

.7 

0 

0 

0 

0 

0 

0 

-.000417 

.002500 

.004060 

.032325 

.8 

0 

0 

0 

0 

■ 0 

0 

0 

0 

-.001523 

.017240 

.9 

0 

0 

a 

0 

0 

0 

0 

0 

-.001077 

.004190 


TABLE III,— VALUES OF THE INTEGRATING .MATRIX [/]' 

(a) Six-Point Solution 


V {2 

0 


.4 


.8 

.g 

0 

0 

0 

0 

0 ^ 

0 

0 

.2 

.08333 

.13333 

-. 01667 

0 

0 

0 

.4 

.06667 

.26607’ 

.06667 

0 

0 

0 

.6 

.06667 

.26007 

.15000 

. 1333r 

-.01067 

0 

.8 

.06607 

.26607 

.13333 

.26667 

.00067 

0 

.0 

.06867 

.26607 ' 

.13333 

.26067 

.10833 

.06607 


(b) Ten-Point Solution 


6'/2 

0 

,1' 

.2 

■ .8 

.4 

.5 _ 

.6 

.7 

A 

.9 

0 

t 0 

0 

0 

0. 

0 

0 

0 

0- 

0 

0 

.1 

.04167 

.00067 

—.00833 

O' 

b 

0 * 

0 

' 0 

0 

0 

.2 

.03333 ' 

.13333 

' .03338 

0 

0 

0 

0 

0 

0 

0 

.3 

.03333 

.13333 

.07500 

.06607 

-.00833 

0 

0 

0 

0 

0 

.4 

.08333 

.13333 

.06067 

.13333 

.03333 

0 ^ 

0 

0 

0 

0 

.5 

.03333 

-.13333 

.06607 

.13333 

.07500 

.00667 ^ 

-.00833 

0. 

0 

Q 

.0 

.03383 . 

.13333 

.00607 

.13333 

■ .06667 

.13333 “ 

.03333 

0 

0 

0 

. 7 

.03333 

.13333 

.00667 

: 13333 

.06607 

. 13333 

.07500 

.06067 

-.00833 

0 

.8 

.03333 

.13383 

.00067 

.13333 

. . 00667 

. 18333 

.06067 

.13333 

.03333 

0 

,0 

.03333 . 

.13333 

.00667 

.13333 

.06607 

. 13833 ^ ; 

.00667 

.13333 

.07500 

.00007 
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TABLE IV.— FORir FOR COMPUTATION OF AUXILIART AEROELASTIC AND ALLERON-REVERSAL MATRICES 

(a) PABAilETERS PbETETBNT TO La.TBEAL^CO’TTBOL CiLCrt-ATIOKS 
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TABLE V.— FORM FOR SOLUTION OF AEROELASTIC EtJUATION FOR LATERAL CONTROL 

(a) Rbversai. (b) Control Power 
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TABLE VI— COMPUTATION OF AUXILIARY AEEOELASTIC AND AILERON-REVERSAL MATRICES FOR THE EXAMPLE 

WING AT SUBSONIC SPEEDS 


(a) Paeametehs Pebtinext to Lateeal-Cokthol Calcceatioks 
6^2=°^ ^2 = 1D00 


£Ej= 0.547 




0-250 -0.2020 -1.0000 


.250 -.1960 -1.0000 


.285 I .458 .Qtll .0563 


1.000 .458 I .0138 | .070L 


LOGO I .458 


LOOO I .453 j .0172 .0906 


(b) CAI.CFLATION: OF THE ALiiLIABT AeHOELASTIC ilATEIX 




0 


imm. 

0 

0 

.00414 

-ooen 

.02541 

.08337 

.03342 

.00774 

-.021S1 

.04563 

.21248 

.00966 

.00774 

-.02181 

.04041 1 

.2S6QI 

.17065 

.00774 

-.02181 [ 

.03443 

.31802 

.20160 

.00774 

-.02181 j 

.03443 

.31802 

.19534 


(c) Calculatiox of the Aileros^ Reversal Matrix . 




0.04325 

Ql 11456 

0.04540 



■ 


(f-UpiJ 

ff=0.2S07 

Lil 

f 


; l"'J 1 

1 



0lI925 0.172S 0.4570 0.1815 0.2158 


¥kl i®' 


0 


0 


.0457 

.1211 

.1725 

.4570 

.1725 

.4570 

.1728 

.4870 



.4 

.6 

0 

0 

.02543 

.02357 

.09385 

.IQ7U 

.15508 

.26327 

.18440 

.38811 

.10101 

.40863 


.00372 


.00665 


-.00732 


.00400 


.04635 





























































































































































